Cotunneling thermopower of single electron transistors 
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We study the thermopower of a quantum dot weakly coupled to two reservoirs by tunnel junctions. 
At low temperatures the transport through the dot is suppressed by charging effects (Coulomb 
blockade). As a result the thermopower shows an oscillatory dependence on the gate voltage. 
We study this dependence in the limit of low temperatures where the transport through the dot 
is dominated by the processes of inelastic cotunneling. We also obtain a crossover formula for 
intermediate temperatures which connects our cotunneling results to the known sawtooth behavior 
in the sequential tunneling regime. As the temperature is lowered, the amplitude of thermopower 
oscillations increases, and their shape changes qualitatively. 

PACS numbers: 73.23.Hk, 73.50.Lw, 72.15.Jf 



CZ2 



-a 
c 

o 
o 



> 
in 

(N 



i 

C 

o 
o 



X 



I. INTRODUCTION 

In the last few years a number of experiments have 
been performed in order to investigate the transport of 
electrons through small conductors, such as-metallic par- 
ticles and semiconductor quantum dots.uu One of the 
most commonly studied types of devices is the single- 
electron transistor which consists of a quantum dot (or 
a small metal particle) connected to two leads by tunnel 
junctions. The particle is usually also capacitively cou- 
pled to an additional gate electrode, Fig. |l|. The trans- 
port of electrons through the quantum dot is strongly 
affected by charging effects. Indeed, when an electron 
tunnels into the dot from a lead, the electrostatic energy 



of the system increases by ~ Ec 



'/2C, where e is the 



elementary charge, and C is the capacitance of the dot. 
In a typical experiment the temperature is low, T Ec, 
and the tunneling is strongly suppressed as only a very 
small fraction of electrons have energy of the order of Ec 
necessary for the tunneling to occur. This phenomenon 
is commonly referred to as the Coulomb blockade. In 
a single electron transistor, the charging energy can be 
controlled by the gate voltage V g . For instance, by ap- 
plying a positive voltage to the gate, one can lower the 
increase in electrostatic energy caused by adding an ex- 
tra electron to the dot. As a result at certain values of 
V g the electrostatic gap vanishes, and the transport of 
electrons through the system is strongly enhanced. This 
leads to a sequence of periodic peaks in the conductance 
of the single electron transistor as a function of V g , which 
is often observed in the experiments.! 2 ! The positions and 
the shape of the peaks are in good affpeement with the 
so-called sequential tunneling theorjocl of conductance 
which accounts for real processes of electron tunneling 
between the dot and the leads. 

In a number of more recent experiments the ther- 
mopower nSi-jOf a single electron transistor has been 
studied.B'BLH3 Similarly to the case of conductance, the 
thermopower shows periodic oscillations as a function of 
the gate voltage. The theory of the thermopower in a 
single electron transistor in the framework of the sequen- 



tial tunneling approach was developed by Beenakker and 
Staringu It was confirmed by the experimental data of 
Ref. ||, where the sawtooth shape of the thermopower os- 
cillations predicted in Ref. j^jvas observed. On the other 
hand, in a recent experirncntfl the observed thermopower 
oscillations had a shape different from the sawtooth, and 
also the amplitude of the oscillations was far smaller than 
the one predicted in Ref. [| 

The most significant difference between the experi- 
ments of Ref. H and Ref. fj] was that the temperature in 
the latter work was very small, T ~ 0.006 e 2 /C, in com- 
parison with the estimated temperature T ~ 0.065 e 2 /C 
in Ref. |j| The authors of Ref, [?] attributed the deviation 
of their data from the theoryu to the fact that the theory 
neglected the effects of virtual tunneling (cotunneling) 
of electrons through the dot. Indeed, it is known that in. 
the case of very low temperatures cotunneling processesEll 
give dominant contribution to the conductance of single 
electron transistors. One can therefore expect that this 




FIG. 1: Thermopower measurement in a single electron tran- 
sistor. The quantum dot is coupled to the two leads by tunnel 
barriers. The electrostatics of the system is controlled by the 
voltage V g applied to a gate coupled capacitively to the quan- 
tum dot. The left and right leads are held at temperatures 
T + AT and T, respectively. As a result a small voltage V 
is generated between the leads. The resulting thermopower 
S = —V/AT is measured as a function of V g . 
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mechanism will result in a different behavior of the ther- 
mopower at low temperatures. 

In this paper we develop the theory of the thermopower 
of single electron transistors in the regime of inelastic co- 
tunneling. This mechanism dominates the low temper- 
ature electron transport in the case of relatively large 
dots, where the effects of finite quantum level spacing 
can be neglected. We find the thermopower oscillations 
of the shape qualitatively similar to that observed in 
the experiment ,U and the amplitude of the oscillations 
is significantly lower than the result of the sequential 
tunneling theory.B We also study the crossover from the 
cotunneling behavior of the thermopower to the saw- 
tooth regime of sequential tunneling which occurs as one 
raises the temperature T of the system above a certain 
crossover temperature T c . The latter is found to be 
T c ~ E c /He 2 /h(Gi + G r )], where G t and G r are the 
conductances of the tunneling junctions connecting the 
quantum dot to the left and right leads. Throughout the 
paper Gi and G r are assumed to be small, Gi,G r -C e 2 /h. 
The opposite regime of strong tunneling, more relevant 
for the conditions of the experiment Ja was recently ad- 
dressed in Ref. [l0| 

In the next section we introduce the theoretical model 
of a single electron transistor used in this paper and dis- 
cuss the relevant mechanisms of electronic transport in 
the device. In Sec. Ill we review the known results for 



the thermopower S in the regime of sequential tunnel- 
ing and obtain S in the regime of inelastic cotunneling. 
In Sec. IV these results are unified in a single formula 
that correctly describes the crossover between the two 
regimes. 



II. MECHANISMS OF TRANSPORT 

To describe the single electron transistor Fig. [j], we 
introduce the Hamiltonian in the form H = Hq + V\ 
where 

Hp = £k%a>k + ^2 £p a t a p + ^2 ^i a \ a i 



2C 



4>Q, 



v = X! (*fc.pi4 a pi + t *k,pi a li 
k,pi 

+ X] (tp2,q a p 2 a q + tp 2 ,q a q a : 
P2,q 

Q = - e ^2 [ a l a P ~ 0(-£p)] • 



(1) 

(2) 
(3) 



Here a^, a p , and a q are the annihilation operators for 
the electrons in the left lead, quantum dot, and the right 
lead, respectively; the electron energies £ p , and £ g 
are measured from the Fermi level; the tunneling in and 
out of the dot is described by matrix elements tk, Pl and 
tp2,q- Operator Q represents the charge of the quantum 



dot, e is the elementary charge, and the potential <f> is 
proportional to the gate voltage V g . 

Our goal is to find the thermopower S of the system 
as a function of the potential 4> for a given temperature. 
The thermopower is defined in terms of the voltage V 
generated across the single electron transistor when the 
temperatures of the left and right leads T\ and T r differ 
by AT <C Ti, T r and no current I through the system is 
allowed: 



S = 



lim — — 

AT^O AT 



(4) 



7=0 



In the linear response regime the current can be presented 
as 



I = G T AT + GV, 



(5) 



where G is the conductance of the system, and the ki- 
netic coefficient Gt describes the current response to the 
temperature difference. The condition I — then results 
in the following expression for the thermopower: 



Gt 
~G' 



(6) 



Thus, one can find the thermopower S by calculating the 
kinetic coefficients Gt and G. 

In this paper we concentrate on the case of relatively 
large quantum dots, where the effects of finite quantum 
level spacing 5 in the dot can be ignored. In the limit 
6^0 the transport of electrons in single electron transis- 
tors can be accomplished via either sequential tunneling 
or inelastic cotunneling mechanisms. The respective con- 
tributions to the linear conductance G are well known. 
The conductance as a function of the gate voltage tj> shows 
periodic peaks centered at points 



>N 



c 



N 



(7) 



where the charging energies of the dot with N and N + l 
additional electrons are equal. The sequential tunneling 
theoryuo accounts for the real events of electron tunnel- 
ing between the leads and the dot. When the (JV + l)-st 
electron tunnels into the dot, the charging energy changes 
by e(4>N — 4>)- At low temperatures T <C e\4>N — <p\, the 
density of electrons with energy sufficient to charge the 
dot is exponentially small, resulting in the conductance 
peaks with exponential tails: 



QSq = 



G]G r 



e(4> - </>jv)/T 



2{Gi + G r ) sinh [e(0-0jv)/T] 



(8) 



On the -other hand, the inelastic cotunneling 
mechanismllll accounts for the second-order tunnel- 
ing processes when, e.g., an electron tunnels from the 
left lead into the dot, and then another electron tunnels 
from the dot to the right lead. The initial and final state 
of such a process have the same charge in the dot. The 
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Coulomb blockade only affects the energy of the virtual 
state, resulting in only a power-law suppression of the 
conductance at low temperatures: 



III. THE THERMOPOWER IN THE REGIMES 
OF SEQUENTIAL TUNNELING AND 
COTUNNELING 



it n 

GCO r~i 

O 6 



T 



[e{(j> - <j) N )Y 



(9) 



The above expression formally diverges at 4> = 0iVi 
because the calculation in Ref. [ll] neglected the con- 
tribution of the quasiparticle energies to the energy of 
the virtual state. To estimate the cotunneling contri- 
bution at the center of a peak, i.e., at 4> = 4>n, one 
can replace the charging energy difference e(</> — c/)n) 
in the denominator of Eq. (ph by the temperature T. 
This results in G co ~ (h/e*)GiG r . The contribution 
(|J) of the sequential tunneling mechanism at the peak 
is G sq ~ GiG r j (Gi + G r ). Thus, near the centers of the 
peaks the conductance is dominated by sequential tun- 
neling. 

On the other hand, between the peaks the sequen- 
tial tunneling contribution (^) decays exponentially at 
T — > 0, as opposed to the relatively slow dependence 
G co tx T 2 of the cotunneling conductance (||). Therefore, 
at low enough temperatures the cotunneling mechanism 
will dominate the conduction in the valleys between the 
peaks. By comparing the contributions (||) and (||) in 
the middle of a valley, i.e., for <j> — 4>n = e/2C, we find 
that the sequential tunneling dominates the conductance 
at any gate voltage only at temperatures T > T c , where 



Ec 



ln[e 2 /n(Gi+G r )Y 



(10) 



At lower temperatures, T < T c , the conductance in the 
regions of width 



T . 



2C T c 



In 



HGl + G r ) 



(11) 



around each peak is determined by the sequential tun- 
neling processes, Eq. (g), whereas outside those regions 
the cotunneling processes dominate, Eq. (g). Note, that 
in our case of weak tunneling, Gi,G r <C e 2 /h, the re- 
gion ([□]) is wider than the thermal width of the peak 
T/e. Thus, the cotunneling mechanism becomes impor- 
tant only away from the peaks, where the conductance is 
much smaller than its peak value. 

The change in the transport mechanism from sequen- 
tial tunneling to cotunneling is more dramatic in the 
case of the thermopower. Unlike the conductance, the 
thermopower in the sequential tunneling regime does not 
have the form of sharp peaks near (f> = <j>N', in fact it 
reaches its maximum values near the centers of the val- 
leys between the peaks of G(<f>), Ref. |[ At T < T c the 
transport in those regions is strongly affected by the co- 
tunneling mechanism, which leads to qualitative changes 
in both the amplitude and shape of the Coulomb block- 
ade oscillations of the thermopower. 



The thermopower S of a system is a direct measure 
of the average energy that the electrons carry during 
the tunpeling processes. This follows from an Onsager 
rclationlla between the Peltier coefficient IT = — (£)/e and 
the thermopower: 



Q _n_ (0 
b ~T~ ir' 



(12) 



The two transport mechanisms in single electron tran- 
sistors, the sequential tunneling and the cotunneling, in- 
volve electrons with typical energies £ ~ Ec and £ ~ T, 
respectively. Thus the thermopower ( |l2| ) is strongly 
affected by the crossover between the two regimes at 
T - T r . 



A. Sequential tunneling regime 

In this section we review the results for the ther- 
mopower in the sequential tunneling regime found by 
Beenakker and StaringJj The sequential tunneling cur- 
rent is determined by the following two elementary real 
transition processes: (a) an electron tunnels between the 
left lead and the dot; (b) an electron tunnels between 
the dot and the right lead. Both processes are of the first 
order in the perturbation parameter Gi ir /(e 2 /h). The 
current in the stationary state can then be derived by 
means of a kinetic equation which involves the probabil- 
ities for the system being in a certain charge state and 
the tunneling rates obtained by Fermi's Golden Rule. 

In the linear-response regime the current can be due 
to either a bias voltage or a temperature difference be- 
tween the two leads. The conductance G and the kinetic 
coefficient Gt are found in Refs. || and ^, respectively: 



s-isq 

\JTrj. 



GlG r 



Gi + G, 



■E< 0) / 



/ En — Eff—i 



N 

1 GlG r 

YeGi+G r 



(0) 



E< 0) / 



N 



T 

En 



(13) 



En -i 



T 

En — En -i 
T ' 



(14) 



Here En = EcN 2 — Ne<j) is the electrostatic energy of 
the dot containing N electrons, f(x) = x/(l — e~ x ), and 
stands for the equilibrium probability distribution 
of the dot charge, W ( ° ] = e - EN ' T / J^n ^ En/t . 

At low temperatures T <C Ec at most two charge 
states contribute significantly to the sums in (|l3|) and 
@. One can then neglect the exponentially small con- 
tributions of the other charge states and obtain Eq. (||) 
for the conductance G as well as an analogous expression 



4 



for Gt, 



GlG r 



2 /T 2 



4(Gi + G r )wah[e{(f>-<f> N )/T\' 



(15) 



Using Eq. (g) one then obtains the low temperature limit 
of the thermopower 



gsq = 



■>N 



2T 



(16) 



Here, cj>N defines the position (Q) of the conductance peak 
closest to 4>. This is the sawtooth like behavior described 
in Ref. ^. The amplitude of the Coulomb blockade oscil- 
lations of the thermopower is given by S^ ax = e/(4CT). 

This result can be understood in terms of the average 
energy of tunneling electrons. Exactly in the middle be- 
tween two conductance peaks (e.g., 4> — [4>n-i + <A/v]/2) 
the same amount Eq of energy is required to either add 
or remove an electron from the dot. Therefore, the two 
processes involving electrons with energies £ « Eq and 
£ « —Eq contribute equally to the electronic transport, 
and the average energy of tunneling electrons (£) van- 
ishes. However, if the gate potential <f) is slightly changed 
by ~ T/e towards one of the conductance peaks, one of 
the two charge states is exponentially suppressed. In this 
case the average energy is found to be (£) as ±Ec- This 
means that the thermopower exhibits a sharp jump be- 
tween two conductance peaks. If the gate voltage is tuned 
further towards </>jv, the average tunneling electron en- 
ergy is given by (£) ~ v,\ = E N+1 - E N = e((j) N - <f>) and 
thus changes linearly with <j). Using this qualitative con- 
siderations the result (16) is reproduced except for an ad- 



ditional numerical factor of 1/2. This factor 1/2 is due to 
the fact that not only electrons with energy u\ above the 
Fermi level are involved in the transport. The probabil- 
ity for an electron with energy £ to tunnel into the dot is 
proportional to the average occupation number ~ e~^/ T 
in the lead and to the density of holes ~ e~^ Ul ~^'' T in 
the dot. Therefore, all electrons with energies less than 
u\ above the Fermi level have the same probability to 
tunnel into the dot. For that reason, one obtains for the 
average tunneling electron energy (£) = u±/2 which then 
reproduces the result (BT 



B. Cotunneling regime 

As we discussed in Sec. [n[ at temperatures below the 
crossover temperature (|i(i| ) cotunneling is the main mech- 
anism of transport in a single-electron transistor. A co- 
tunneling process consists of two steps. First, an electron 
close to the Fermi level in one lead tunnels into the dot, 
thus changing the energy by ~ u\. Since the energy is 
not conserved during this process, the system is in a vir- 
tual intermediate state. A second tunneling process has 
to follow which involves an electron tunneling from the 
dot to the other lead and thus restoring the conservation 
of energy. Since the charge of the dot before and after a 




Pi 



FIG. 2: The two types of inelastic cotunneling processes 
transferring electrons from the left to the right lead. The first 
process consists of the following two steps: an electron tunnels 
from a state k in the left lead to a state pi in the dot, and 
then an electron from a state pi in the dot tunnels to state 
q in the right lead. The second process transfers electrons in 
the opposite order: it starts with an electron tunneling from 
the dot to the right lead, and finishes with another electron 
tunneling from the left lead into the dot. The energies of the 
virtual states in the two types of processes are affected by the 
electrostatic energies ui and u_i required to either increase 
or decrease the charge of the dot by that of one electron. The 
dash-dotted lines show the positions of the Fermi level in the 
leads and the dot. 



cotunneling process is the same, the electronic transport 
is primarily due to electrons with energies £ ~ T. There- 
fore, unlike sequential tunneling, the cotunneling contri- 
bution to the current is not exponentially suppressed at 
T <C Ec , and it dominates the transport in the valleys of 
the Coulomb blockade. In the following only inelastic co- 
tunneling processes which involve two different electrons 
are considered, since the elastic ones are suppressed at 
small level spacing in the dot, S <C T 2 /E c , Ref. |l[ 

The cotunneling contribution to the current is of the 
second order in the small parameter Gi :T / '(e 2 /Ti) because 
two coherent tunneling processes have to take place. For 
a given number N of electrons in the dot there are two 
types of second order processes which can transfer elec- 
tron from one lead to the other, see Fig. 0. The current 
then reads 



2tt 



^2 [ n k0- - n pi )n P2 (l - n q ) 



-n q (l -n P2 )n Pl (l - n k )]S(£, k - £ Pl + £ P2 - f 9 ) 



5 



P2,g 



pi b P2,q 



Ul 



6 - Cpi + u -l 



(17) 



where n^, n p 



and n„ are the occupation numbers 



of the respective states; u±\ = Em±i — -Bjv, and N is 
determined by minimizing the electrostatic energy of the 
dot for a given potential <fi. Note that the first product of 
the occupation numbers corresponds to the current from 
the left to the right lead, whereas the second one accounts 
for the current in the opposite direction. 

If the bias voltage V is zero, and all three electrodes 
are at the same temperature T, the current from the left 
lead to the right lead cancels with the one flowing into 
the opposite direction. This can most easily be seen by 
exchanging k <-» q and p\ <-> p2 in one of the products 
of occupation numbers which is possible because of the 
apparent left-right symmetry of the system. 

The cotunneling contribution G co to the conductance 
can be derived by linearizing Eq. (|l7]) with respect to the 
bias voltage V. After replacing the sums by integrals and 
the occupation numbers by Fermi functions, the conduc- 
tance can be presented as a single integral, 



TlGlG r 

8ttc 2 T 



e 



sinh 2 



1 



d£. (18) 



This expression formally diverges at £ = ±u±i; w e wil l 



discuss the proper regularization procedure in Sec. IV A 



At low temperatures T <C the contributions of 

the regions near the singularities are exponentially small, 
and one can evaluate G co by neglecting the quasiparticle 
energies £ ~ T in .the denominators. Then the conduc- 
tance G co is foundlld to be 



G c ° = ^G l G r T 2 



1 



(19) 



If the potential (j> is significantly closer to the peak at <pN 
than to 4>n—i the term 1/u—i can be neglected compared 
to 1/tti, and Eq. (^|) is reproduced. 

In order to calculate the kinetic coefficient Gt we as- 
sume a slightly higher temperature T + AT in the left 
lead compared to the dot and the right lead. Then the 
current ( |l7|) is linearized with respect to the small tem- 
perature difference AT <C T. After replacing the sums 
by integrals and occupation numbers by Fermi functions, 
the following expression for Gt is obtained: 



Geo 
T — 



TlGlGr 

27re 3 
x T d£ k 



J d£k d£ Pl d£ P2 d£ q 5(£ k - f Pl 



£p2 £<? 



1 



1 



{[1 - n(£ pl )]n(£ P2 )[l - nfo)] + n(&J[l - ra(£ Pa )]n(£,)} . 

I 



(20) 



The same expression is obtained in a more careful treat- 
ment where the temperature of the dot is not necessarily 
equal to the temperature of the right lead. Three of the 
four integrals can be calculated exactly, resulting in the 
following expression for Gt 



_ HGiG r 1 
° T " Wire 3 T 2 



sinh 2 ^ 



1 



1 



(21) 

As it was the case for the conductance G, the main con- 
tribution to Gt comes from energies £ of the order of 
T <C Ux,U-%. However, setting those terms in the de- 
nominators to zero yields a vanishing Gt as the inte- 
grand in (21) becomes an odd function in £. The first 
non- vanishing contribution to Gt is obtained by expand- 
ing the fractions in the integral up to first order in £,/u\. 
The final result for G§? can be written as 



4tt 3 h 



GCO *" f~i rp3 I 
rp — — Lz / Lx r 1 J 



(22) 

Using the expression (0) for the thermopower we find 



from (H) and © 



5 e 2 \(f> - cj) N (j>- 



(23) 



This expression corresponds to the potential range be- 
tween two conductance peaks 4>n~i < 4> < 4>n with 
the exception of the values of <fi very close to the peaks, 
namely (4> N - 0), (0 - 4>n-i) > T/e. 

The result (|2^) can be understood in terms of the av- 
erage energy of tunneling electrons, Eq. (jl|). The trans- 
port in the cotunneling regime is mainly due to the elec- 
trons with energies £ ~ ±T. In the low temperature 
limit T <C Mi the tunneling probability from a state with 
energy £ into an intermediate state with energy £ p is pro- 
portional to 



u>(£) oc 



1 



1 



111 



(24) 



This expression clearly shows that the tunneling prob- 
ability w(£) is enhanced for electrons above the Fermi 
level £ > 0, and therefore the average energy (£) will not 
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vanish. In fact, with a typical energy £ ~ T for electrons 
involved in the tunneling and the probability ( |24| ) we find 
for the average energy (£) ~ T 2 /ui. Using the relation 
(12), the cotunneling contribution to the thermopower 
(23) is reproduced in the correct order of magnitude. 



IV. THE THERMOPOWER AT ARBITRARY 
GATE VOLTAGE 

In the previous section expressions for the ther- 
mopower in both the sequential tunneling and cotunnel- 
ing regime were presented. The sequential tunneling re- 
sult is given by Eq. ( |l6| ) . It dominates the thermopower 
for gate potentials close to the positions of the conduc- 
tance peaks 4>m- On the other hand, the cotunneling re- 
sult (23) gives the correct description of the thermopower 
between the conductance peaks 4>n-i < 4> < 4>N if the 
temperature T is below the crossover temperature T c . In 
this section we find the combined contributions of both 
transport mechanisms and obtain an expression for the 
thermopower valid at any gate voltage. 

The current through the quantum dot is the sum of 
the two contributions, sequential tunneling and cotun- 
neling. Therefore using Eq. (||) and (||) we find for the 
thermopower 



S 



G s i + G c 



(25) 



The well knowrOJ sequential tunneling contributions G sq 
and are valid for the entire potential range. How- 
ever, the cotunneling results ([l9]) and ((2^) are only 
valid away from the centers of the conductance peaks, 
i.e. (<f>N — <j>), {4> — (j>N-i) ^> T/e, and diverge close to 
the peaks. In the following we will show how to cor- 
rectly regularize these divergences. First we discuss the 
thermopower (|2^) in the limit of very low temperatures 
T <C T c . Then we present the more general result valid 
for higher temperatures. 



The low temperature thermopower in the 
vicinity of a conductance peak 



part r to the energy of the intermediate virtual states in 
Eq. @, 



k,pi,P2,q 



n k (l - n pi )n P2 (l - n q ) 



n q (l -n P2 )n Pl (l -n k )]S(^ k -£ Pl + £ P2 - £,) 

2 



ik - £ P i -ui — iT 



(26) 



In the limit r — > this expression can be divided into 
a large part ~ 1/r and a smaller one that is independent 
of r using the following procedure: 



dE 



f(E) 

e 2 + t 2 



dE 



/(0) 



E 2 



dE 



f(E)-f(0) 



r 2 J e 2 + y 2 

f(E) - /(0) 



dE- 



E 2 



(27) 



where the last integral is to be understood as the prin- 
cipal value. As we show in the Appendix |X|, the hrst 
term corresponds to the sequential tunneling. In partic- 
ular, the results (||) and ( |l5| ) are restored in th e li near 
response regime using the scheme ( |2~7j ) and Eq. (A6) for 
r. The second term represents the correct non-divergent 
expression for the cotunneling. We apply the scheme j27| ) 
by rewriting the cotunneling conductance ( |l8| ) and the 
kinetic coefficient (^l|) according to the following rule: 



dE 



1(E) 
E 2 



dE 



f(E)-f(0) 
E 2 



dE m+f{-E)-2m) (28) 



This regularization procedure coincides with the one dis- 
cussed in Ref. [l4] for the case of multichannel tunneling 
junctions. The result for G co and G§? then reads 



h 



Aire 2 

h 



' 47re 3 



GiG r T 

GlG r Tt 



2T 

e(4> - 4>n) 
2T 



(29) 
(30) 



where the functions T and Tt are defined as 



In the regime of very low temperatures T <sC.T c there 
are at most two charge states, e.g., N and N + 1 elec- 
trons in the dot, that contribute significantly to the to- 
tal current while the contributions of all other states 
are exponentially suppressed. The crossover between the 
sequential tunneling and the cotunneling regimes takes 
place at a gate potential close to the conductance peak 
\4> — 4>n\ ~ A0, where A</> is given by Eq. (O). In order 
to find the behavior of the thermopower (P5|)in this po- 
tential range we can neglect terms ~ occurring in 
the current (|i7|). 

To regularize the singularities at — £ Pl = u\ in the 
expression (|l^) for the cotunneling current we follow the 
approach presented in Ref. O and add a small imaginary 



dz 



(1 



sinh [x(l + z) 



(1-z) 2 



sinh [.t(1 — z)\ sinh [a 



xx 



dz 



(l + z) 3 



sinh 2 [a; (1 + , 
sinh 2 [a;(l — z)\ 



sinh 2 fxl 



,(31) 



.(32) 



The final result for the thermopower is then obtained 
by substituting the sequential tunneling results (0) and 
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A<|>/(T/e) 



FIG. 3: The thermopower S of a single electron transistor in 
the vicinity of a conductance peak. The distance from the 
peak A<f> is assumed to be small compared to the distance 
to the next peak, A<f> < e/C. The solid curve is calculated 
for Gi + G ,- = 10 -3 x by substituting the expressions 

(|), (|l|), (§|), and @ into Eq. (§|). The asymptotics of 
cotunneling and sequential tunneling are shown by dashed 
and dash-dotted lines, respectively. 



eA(t> 2 



In 



- (in- 

9 V 9 



(34) 



The thermopower at A(f> < A(f>i is given by the sequen- 
tial tunneling result S ~ Acf>/2T. In the narrow range 
A0i < A(f> < A(f>2 the conductance is already dominated 
by the cotunneling, while the main contribution to Gt 
in the numerator of (|2^) is still due to the sequential 
tunneling. Therefore we find a steeply descending ther- 
mopower S oc exp[— eA<p/T] in that region. Eventually, 
at <]) > 02 both the numerator and the denominator in 
Eq. ( p5| ) are governed by the cotunneling contribution, 
and S ~ T / (eA(f). Therefore the maximum of the ther- 
mopower is reached at A<fi ~ A0i and can be obtained 
using the sequential tunneling expression (plf): 



S -ilni 

2e g 



(35) 



To find the value of the thermopower after the exponen- 
tial fall-off, one can substitute A</>2 into the cotunneling 



expression (fc 

Smax - 



which gives S ~ (47r 2 /5)e 
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( |l5| ) and the regularized cotunneling results (g9|) and ( J30| ) 
into Eq. (|2^) . The dependence of the thermopower on the 
gate voltage obtained from those expressions is illustrated 
in Fig. |. 

As can be seen in Fig. y, the thermopower is domi- 
nated by the sequential tunneling contribution for small 
Acf> — 4> — 4>n. After it reaches its maximum value S max , 
the thermopower falls off sharply to merge with the co- 
tunneling result (^3|) for large A(f>. To understand this 
behavior one should notice that the crossover between 
the sequential tunneling and the cotunneling occurs at 
different values A<fii and A</>2 for the conductance G and 
for Gt- The two crossover values of the gate voltage 
can be estimated by setting G sq = G co and G^ = Gff. 
At small conductance g = h(Gi + G r ) / (2-Ke 2 ) <C I the 
Eqs. (|), ©, © and (H) result in 



B. The thermopower at arbitrary temperatures 



eA^i 
T 



In 



- (in- 

5 V 9 



(33) 



In Sec. [VA we studied the thermopower at very low 
temperatures T <C T c . In this regime the thermopower 
in the valleys between the conductance peaks centered 
at <f> = 4>n is described by the simple cotunneling result 
( p3| ) ; the sequential tunneling contributes only in narrow 
regions around <fi = (f>j\f. Another interesting regime is 
that of the temperatures of the order of T c , where one can 
explore the crossover from the sawtooth behavior ([ll]) to 
the correct low-temperature limit. 

The crossover takes places in the valleys between the 
peaks, where more than two charge states give compara- 
ble contributions to the current through the device. Thus 
we present the total cotunneling current as a weighted 
sum of the contributions of states with a different num- 
ber N of electrons in the dot, I = ^2 • 
is given by 



'Iff, where Iff 



2tt 



Iff = — — e [ n k{l-npJnp 2 (l-n q )-n q (l-n P2 )n pl (l~n k )}5(t k -£, Pl +£, P2 -£ q ) 



Pi u P2,q 



Pi b P2,q 



£fc — Cpi — + En £fc — £ Pl + E^-i — En 



(36) 



Note that one has to keep terms ~ f/it-i = 1/(En-i — Em) if the expression is to be valid not only close to the 
peaks at 4>n but also in the valleys between the peaks. Similarly to the discussion in the previous section, we linearize 
the current with respect to a small bias voltage V and to a small temperature difference AT, yielding G and Gt- 
The cotunneling contributions have then to be regularized applying the rule (Eq) to each of the divergent terms. The 
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result is 

QCO 



h GiG r 



e 2 47r 



GCO 
T — 



e 3 Air 



E 

— c 
o 

E 



oo 



r(0h 



Em — E 



N- 



where the functions J-* and J-j, arc defined by 



F*(x) = x\x\ 



dz 



oo 

T* T {x) = x 2 \x\ J ^ 



2T 

En — Eff-i 
2T 



(1 + zf 



{1-zf 



z \sinlr[a;(l + z)} smh 2 [x(l - z)} 



(l + z) 3 



(l-zf 



sinh [x(l + z)] sinh [x(l — z) 
I 



2T 

En — En-i 
2T 



(37) 
(38) 



(39) 
(40) 



The nondivergent cotunneling contributions G co and 
together with the well known sequential tunneling re- 
sults (Q and (Q) for G s « and Gy 9 then give the ther- 
mopower according to Eq. (p5|). This expression for the 
thermopower S 1 is the final result of the calculation. It is 
valid at any gate voltage 4>. As the temperature changes, 
this result shows the crossover from the sawtooth shape 
of the Coulomb blockade oscillations of the thermopower 
at T 3> T c to the low-temperature behavior discussed in 
Sec. [TV A| at T < T c , see Fig. |. The low-temperature 
curve of S{4>) shown by dotted line in Fig. [| is qual- 
itatively similar to the experimentally measured ther- 
mopower of Ref. 0. 



— T/E c 


= 0.1 


~~ T/E c 


= 0.05 


T/E c 


= 0.025 



S/e 




FIG. 4: The thermopower of a single electron transistor 
at different temperatures. Below the crossover temperature 
T c the thermopower in the valleys of Coulomb blockade is 
suppressed, and the shape of the oscillations changes from 
the sawtooth (solid line) typical of sequential tunneling to the 
low-temperature behavior shown by dotted line. The curves 



V. SUMMARY 



We have studied the thermopower of a single electron 
transistor based on a quantum dot weakly coupled to 
two leads, Fig. |l|. The transport through the dot is gov- 
erned by two different electronic transport mechanisms: 
the sequential tunneling and the cotunneling. At tem- 
peratures above the crossover temperature T c given by 
Eq. (10) the sequential tunneling dominates. The ther- 
mopower in this regime was studied in Ref. |9[ At tem- 
peratures below T c the cotunneling mechanism gives the 
main contribution to the transport in the valleys be- 
tween the Coulomb blockade peaks of conductance. This 
changes the shape of the Coulomb blockade oscillations of 
the conductance dramatically, Fig. ||. We derived a single 
expression for the thermopower at arbitrary temperature 
including both contributions. It is given by Eq. ( |25| ) in 
combination with Eqs. (|||), ((lj), @, and @. As can 
be seen in Fig. ^, the resulting thermopower significantly 
changes its shape from the sawtooth behavior for tem- 
peratures above T c to the low-temperature behavior at 
T<T C . In the regime of low temperatures the amplitude 
of the thermopower becomes temperature independent, 
Eq. (^). Near the maximum of the thermopower the 
transport is determined by the balance of the two mech- 
anisms. The sh ape of the thermopower peaks is studied 
in detail in Sec. IV A and illustrated in Fig. ||. 



were calculated for Gi + G r = 10 x 
Eqs. ©, (0), @, and © into Eq. (| 



by substituting 
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APPENDIX A: SEQUENTIAL TUNNELING 
LIMIT 



In this appendix we show how the sequential tunneling 
contribution I sq can be extracted from the current ( p6| ) 
by taking the Y — > limit and keeping only large terms 
1/r. Using the scheme ( p7f ) the fraction !/[(£& — £ Pl 



Ul) 



T 2 ] is replaced by 7T(5(£fc — £ Pl — ui)/Y. The result 



can be presented as 



jsq 



R : 



N^N+l Rn-\ I. — A 



T}1 nr 



2Y/K 



( A1 ) 

where the tunneling rates R are given by Fermi's Golden 
Rule. For the tunneling between the left lead and the 
dot they are 



R 



k,p 



R 



i 



2tt 



^ \tk, P \ 2 n p (1 - rife) - - wi)j 



k.p 



and i? r can be obtained from the expressions for R l by 
replacing the indices k — ► q. n 
On the other hand, in the sequential tunneling theorycl 
the current is presented as 

I s " = -e(W N R l N ^ N+1 ~W N+1 R l N+1 ^ N ) 

= -e(W N+1 R r N+1 ^ N -W N R r N ^ N+1 ),(A2) 

where Wjy represents the probability of the system be- 
ing in a state with N electrons in the dot. These two 
expressions for the current I aq arise from the condition 



that the current through the left barrier equals the cur- 
rent through th e right barrier. Employing this equality 
the current ( A2 ) can be rewritten as 



jsq 



-eW N 



Ft 1 4- R r 



(A3) 

This is the same relation as Eq. (Al) if we make the 
identification 



N 



(A4) 



This expression coincides with the one proposed in 
Ref. PL 



R 



R 



N + l^N 



R 



N- 



*N+1 + R-N^N+l) 

(A5) 



if one takes into account the balance equation (A2) and 
the condition Wjv + Wjy+i = 1- Our derivation shows 
that replacing the denominator in ( p6| ) with a 5-function 
reproduces the result of the sequential tunneling even if 
eV and/or AT ar e not small compared to the temper- 
ature T. In Sec. IV A we use this result in the linear 



response regime, eV, AT <C T, to conclude that this ap- 
proximation reproduces the results (g) for G sq and ( |l5| ) 
for . In this case the Y is determined by the equilib- 
rium rates and can be calculated by replacing the occu- 
pation numbers with Fermi functions, 



Y = ^{Gi + G r )ui coth 



Ml 



2T 



(A6) 



in agreement with Ref. 
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